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FOREWORD 


This  quarterly  progress  report  covers  the  research 
accomplished  under  contract  Nonr  -30l6(00)  through 
31  August  1962.  The  personnel  associated  with  this 
program  were  M.  V.  Peck,  T.  Houser  and  H.  N.  Chu. 
Dr.  S.  S.  Penner  was  a  consultant. 


ABSTRACT 

The  radius  vs  time  data  taken  at  two  oxygen  flowrates 
were  fitted  by  the  empirical  equation  r  A  +  Bt  . 

The  instantaneous  burning  rate,  i,  is  obtained  by  dif¬ 
ferentiation  of  this  equation.  The  tbersmcouple 
measurements  were  described  by  the  psuedo- steady- 
state  heat  conduction  equation  T  «  Tg  ezp  (-bx/k). 

A  theoretical  development  from  the  von  Karman  fun¬ 
damental  equations  for  aerothermochemistry  did  not 
improve  the  theoretical  model  previously  developed. 
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INTRODUCTION 


This  research  program  in  hybrid  combustion  has  been  concerned  with  the 
investigation  of  the  processes  involved  in  the  reactions  between  solid 
fuels  and  gaseous  oxidizers.  The  purpose  of  these  studies,  in  the 
presence  of  a  flame,  has  been  to  develop  a  suitable  model  for  the  hybrid 
combustion  processes. 

In  previous  reports  (Ref.  1  thru  3)  the  experimental  results  indicated 
that  the  hybrid  combustion  of  plexiglass  and  gaseous  oxygen  was  consist¬ 
ent  with  the  idea  that  the  combustion  rate  is  largely  controlled  by 
diffusion  processes.  A  simplified  theory  of  hybrid  combustion  of  the 
same  order  of  sophistication  as  the  classical  treatment  of  Burke  and 
Schumann  (Sef.  6)  for  the  flame  shape  in  cylindrical  diffusion  flames 
was  presented  for  the  plexiglass-oxygen  system  (Ref.  l).  With  proper 
definition  of  the  diffusion  parameters,  this  theory  may  be  applied  to 
all  types  of  flow  and,  in  particular,  to  the  turbulent  flow  conditions 
examined  experimentally. 

The  purpose  of  the  more  recent  work  has  been  to  define  better  the 
contributions  of  the  experimental  parameters  to  the  observed  regression 
rates  (Ref.  4  and  3)'  In  this  report  period  effort  was  directed  toward 
obtaining  better  defined  functions  of  radius  with  time  at  the  two  en¬ 
trance  Reynolds  numbers  (33*000  and  36,000)  and  from  these  functions 
calculating  instantaneous  regression  rates.  In  addition  the  problem 
of  temperature  profile  in  the  solid  fuel  was  explored  further. 
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EXPERIMENTAL  APPARATUS  AND  PROCEDIEE 

The  basic  apparatus  (.Pig.  .1)  has  been  described  in  previous  reports 
(Ref.  .1  and  5)  The  ’iubular  configuration  is  used  for  the  present 
experimenta.l  program.  The  mixing  and  f.lov  stabilization  tube,  which 
was  a  2-sq  in.  duct,  has  been  replaced  by  a  I-inch-diame ter  tube  to 
prevent  the  turbulence  caused  by  the  abrupt  change  in  diameter  at  the 
entrance.  In  a  typical  run,  a  fuel  cylinder  composed  of  3  by  3  by 
2-inch  blocks  with  1-inch  diameter  concentric  hole.*  aligned  with  the 
^  holes  of  the  facing  and  backing  plates,  is  burned  for  the  specified 

time.  The  axial  length,  the  composition  and  temperature  of  fuel, 
the  time  of  burning,  and  the  initial  diameter  of  the  concentric  hole 
may  be  varied.  Thermocouples  may  be  imbedded  between  adjacent  fuel 
blocks  at  any  desired  radial  distance  to  measure  the  increase  in 
temperature  at  that  position  as  the  burning  surface  approaches. 
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Figure  1.  ^brid  Coabustion  Apparatus 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 

DEPENDENCE  OF  RADIUS  ON  TIME 

To  determine  the  change  in  the  radius  of  the  fuel  surface  with  burning 
time,  runs  have  been  made  at  various  burning  times  with  plexiglass 
(PMMA)  fuel  blocks  and  oxygen  at  two  entrance  Reynolds  numbers  of 
approximately  33|000  and  36.000  The  final  radius  was  calculated 
from  the  weight  burned  for  each  axial  position  at  each  burning  time, 

assuming  symmetrical  burning. 

V 

I  Since  a  line  pressure  of  about  3  or  4  psig  was  observed  at  the  start 

of  combustion  when  Re,.  -  33*000  (which  reduced  to  0  psig  at  10  to 
13  seconds),  a  second  series  of  runs  was  made  at  the  lower  flowrate 
^Reg^  *  36,000)  where  no  initial  line  pressure  was  observed.  Twelve 
runs  were  made  at  this  oxygen  flowrate  using  an  axial  length  of  24  inches. 

Figure  2  is  a  comparison  of  the  two  flowrates,  in  which  a  typical 
plot  of  the  final  radius  vs  burning  time  for  block  3  (axial  position 
of  3  inches)  for  both  Re.  =  33 1 000  and  36,000  is  shown. 

It  was  hoped  that  it  would  be  possible  to  secure  a  definition  of  the 
radius  (r)  as  a  simple  function  of  time  (t)  which  could  be  utilized 
over  the  entire  range  of  time.  The  Research  Multiple  correlation  computer 
program  which  applies  various  statistical  tests  in  attempting  to 
obtain  the  best  fit  of  the  supplied  functions  to  the  experimental  data, 
was  used  to  acquire  the  best  function. 
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Figure  2.  Rediue  of  FMMA  Block  3  vs  Buming  Tiae 
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Among  the  functions  fitted  to  the  r  vs  t  data  was  r  =  A  +  Bt  ,  where 

c  =  0.75  for  t  60  and  c  =  1.0  for  t  ::  60  gave  the  best  fit.  The 

fit  for  a  single  function  over  the  entire  time  range  was  almost  always 

poor  at  the  low  times.  The  curve  fitting  constants  are  summarized  in 

Tables  1  and  2;  the  curves  for  block  3  are  shown  in  Fig  2  Upon 

0  75  -0  25 

differentiation  of  r  =  A  Bt  '  ,  r  was  equal  to  0.75  Bt  ‘  ,  which  is 

infinite  at  time  zero.  Although  this  is  unrealistic  from  a  physical 
standpoint,  it  can  be  compared  to  the  infinite  heat  flux  for  a  semi¬ 
infinite  slab  at  time  zero  (i.e.,  time  zero  is  not  physically  defined 
in  either  case).  Therefore,  the  equation  for  r  does  not  have  to  be 
defined  at  t  =  0  to  make  the  r  function  useful  Other  functions  were  not 
as  good,  as  shown  below. 

2 

The  r  vs  t  data  was  also  fitted  well  by  the  function  t  -  A  -f  Br  4*  Cr 
which  yielded  a  rate  equation  r  =  (B  +  2Cr)'^.  This  equation  is  finite 
at  time  zero.  However,  since  experimentally  the  burning  times  are 
defined  and  then  the  radii  determined,  it  is  statistically  improper  to 
reverse  this  and  define  t  as  the  dependent  variable  and  r  as  the  inde¬ 
pendent  variable. 

2  3 

The  function  r  =  A  +  Bt  +  Ct  +  Dt  appeared  to  be  a  satisfactory  fit 

for  the  r  vs  t  data.  In  general.  A,  B  and  D  were  found  to  be  positive 

and  C  was  negative.  At  high  t  the  calculated  rates  (r)  were  found 

2 

to  increase  because  the  positive  term  3Dt  increases  much  more  rapidly 
than  the  negative  term  2Ct.  This  is  impossible  since  no  inflection 
was  observed  in  the  r  vs  t  plots,  thus  the  fit  was  not  as  good  as 
first  believed. 


TABLE  1 


CURVE  FITTING  CONSTANTS* 


J 

i 

B  X 

100 

55,000 

36 , 000 

55,000 

Block 

1 

1.27 

2.67 

2.35 

2 

1.23 

2.54 

2.68 

3 

1.24 

1.25 

2.36 

2.50 

k 

1.24 

1.28 

2.39 

2.46 

5 

1.24 

1.27 

2.49 

2.62 

6 

1.22 

1.26 

2.60 

2.79 

7 

1.23 

1.26 

2.52 

2.87 

8 

1.23 

1.28 

2.49 

2.83 

9 

1.24 

1.28 

2.42 

2.92 

10 

1.25 

1.30 

2.46 

2.88 

11 

1.26 

1.32 

2.58 

3.15 

12 

1.26 

1.30 

2.99 

2.53  . 

*r=A  +  Bt 


0.75 


for  t  —  60 


i 
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It  should  be  noted  that  the  linear  portions  of  the 
60  seconds  are  the  best  straight  lines  through  the 
tinuations  of  the  curves  below  60  seconds. 


curves  for  t  above 
data  and  are  con- 


INSTANTANEOUS  BURNING  HATES 

The  instantaneous  burning  rates  can  be  tabulated  in  two  ways:  r  vs  z  at 
constant  time  and  f  vs  z  at  constant  r.  It  is  believed  that  the  values 
for  blocks  1  emd  12  exhibit  end  effects  and  so  have  not  been  included  in 
the  tabulations. 


Burning  Rate  as  a  Function  of  Time 

By  differentiation  of  the  empirical  equations  of  r  vs  t,  it  is  possible 
to  obtain  the  instantaneous  rates  for  any  time  and  axial  positions.  Repre¬ 
sentative  values  are  summarized  in  Tables  3  and  4.  The  f  changes  by  a 
factor  of  about  4  at  both  flowrates  over  the  burning  times  used.  At  the 
lower  flowrate,  the  f  does  not  appear  to  change  significantly  with  z, 
whereas  the  higher  flowrate  gives  a  steadily  increasing  f  at  all  times. 
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Burning  Rate  as  a  Function  of  Radius 

At  the  higher  flowrate  the  radius  varies  with  axial  position,  and  there¬ 
fore,  at  any  specified  time,  the  rates  are  not  strictly  comparable  be¬ 
cause  the  Beynolds  number  varies  with  radius  It  is  believed  that  a 
better  basis  of  comparison  of  r  vs  z  would  be  at  equal  radii  rather  than 
at  equal  burning  times.  Table  3  summarizes  tbe  rates  at  three  radii  for 
=  36,000  and  55,000.  The  ratio  of  the  rates  at  the  two. oxygen 
flowrates  (as  given  in  Table  3)  are  summarized  in  Table  6,  and  the  dependence 
on  axial  length  is  observed  for  the  ratio 

At  this  point  the  effects  of  turbulence  on  regression  rates  must  be  con¬ 
sidered  as  indicated  by  these  two  sets  of  data.  From  the  relative  regre- 
sion  rates  at  the  two  oxygen  flowrates,  and  assuming  HgO  and  CO  or  COg 
as  the  products  of  oxidation,  the  relative  volume  rates  of  flow  can  be 
calculated.  It  is  found  that  the  relative  volume  rate  of  flow  changes 
very  slightly  with  z  and,  in  fact,  appears  to  get  smaller;  that  is,  at  the 
entrance  the  ratio  of  oxygen  flowrates  is  about  I.30,  and  at  block  11, 
the  flowrate  ratio  of  oxygen  and  products  is  about  1.42  for  HgO  and 
CO,  and  1.46  for  HgO  and  COg  Although  the  absolute  volume  rate  of  flow 
increases  with  z,  the  relative  flowrate  does  not.  Since  these  flowrates 
are  at  constant  radius,  the  relative  turbulence  must  remain  about 
constant.  The  relative  regression  rates  remain  constant  at  about  1.10 
up  to  10  inches,  and  then  steadily  increase  to  above  1.35,  at  20  inches. 
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TABLE  5 

INSTANTANEOUS  BUBNING  RATES,  r  x  1000  cm/ sec 
Constant  Radii 


Radius,  centimeters  1.30 

1.65 

1.80 

Block  No. 

Rcg^  =  36.000 

o 

55,000 

36,000 

55,000 

36,000 

55,000 

2 

8.7 

9.5 

7.5 

8.2 

6.8 

7.3 

3 

8.0 

8.8 

6.8 

7.5 

6.2 

6.7 

4 

8.1 

8.9 

7.0 

7.5 

6.3 

6.7 

5 

8.5 

9.5 

7.3 

8.1 

6.6 

7.2 

6 

8.8 

10.2 

7.7 

8.7 

6.9 

7.8 

7 

8.6 

10.6 

7.4 

9.0 

6.7 

8.1 

8 

8.4 

10.7 

7.3 

9.0 

6.6 

8.0 

9 

8.2 

11.2 

7.1 

9.4 

6.4 

8.4 

10 

8.5 

11.8 

7.3 

9.8 

6.5 

8.7 

11 

9.2 

13.2 

7.8 

10.8 

7.0 

9.5 
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It  appears  that  the  major  difference  between  the  two  flows  is  that  the 
concentration  of  oxygen  remains  higher  down  the  burner.  This  indicates 
two  possible  explanations  for  the  observed  change  in  relative  rates; 

1.  The  average  gas  temperature  may  remain  higher  at  the  higher  flow 
since  the  combustion  can  be  more  complete,  and/or 

2.  Oxygen  may  reach  the  surface  and  aid  in  the  unzipping  of  the 
polymer;  thus,  as  the  concentration  difference  becomes  greater 
down  the  burner,  the  regression  rate  differences  may  also  increase. 

These  two  possibilities  cannot  be  differentiated  at  this  time  and  both 
may  be  operating. 


THERMOCOUPLE  MEASUREMENTS 

The  results  of  the  thermocouple  measurements  are  summarized  in  Table  7  • 
The  experimental  data  was  found  to  be  fitted  by  the  following  equation: 

T  =  T  exp  (-bx/k) 
s 

2 

k  =  diffusivity  in  solid,  cm  /sec 

X  z  linear  distance  from  thermocouple  to  burning  surface, 
centimeters 

(x  +  r)  =  the  sum  of  radius  and  distance  to  the  surface,  a  constant 
T  K  temperature  at  thermocouple  referenced  to  block  temperature 
T^  s  surface  temperature  referenced  to  block  temperature 
r  =  average  burning  rate ,  cm/ sec 
b  =  [k/(x  +  r)  +  rj 
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^Referenced  to  block  tenperature  of  approziaately  23  C 

‘''Calculated  rates  are  derived  froa  the  slope  of  T  vs  x  curves,  the  observed 
rates  are  the  averages  for  the  run. 

'^'The  axial  distances  of  the  theraocouple  for  these  two  runs  were  18  inches, 
for  all  others  it  was  4  inches. 
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This  equation  was  obtained  using  the  assumptions  that  temperatures  and 
rates  were  independent  of  time  (steady  state  and  constant  regression  rates) 
and  solving  the  moving  boundary  differential  equation  for  a  cylindrical 
surface  (Ref.  5,  Eq.  7): 

St/  St  *  k  S  ^t/S*^  +  [k/(*  +  r)  +  f  ]  S^/S*  =  0 


or 


d^T/dx^  *  -(b/k)dT/dx 

Figure  3  shows  typical  examples  of  the  consistency  of  fitting  this  equa¬ 
tion  to  runs  1431  (39  seconds  burning  time)  and  1432  (41  seconds  burning 
time).  As  indicated  in  Table  7,  the  calculated  and  observed  burning 
rates  are  in  fair  agreement.  This,  coupled  with  the  good  fit  of  the  de¬ 
rived  equation  to  the  temperature-distance  data,  indicates  that  the  in¬ 
correct  assumptions  of  steady  state  and  constant  burning  rate  necessary 
for  the  derivation  of  the  equation  do  not  prohibit  its  use. 

A  discussion  of  the  probable  error  and  the  significance  of  these  measure¬ 
ments  is  in  order.  The  range  of  surface  temperatures  obtained  can  be 
explained  easily  by  the  uncertainty  in  the  measurement  of  x,  which  is 
known  no  better  than  ±0.3  millimeter.  This  is  indicated  by  the  two 
curves  in  Fig.  3>  Since  this  type  of  error  is  due  to  random  uncertainty 
an  average  over  the  14  runs  listed  in  the  table  should  produce  an  ac¬ 
curate  surface  temperature,  assuming  that  a  systematic  error  does  not 
influence  the  data.  It  is  believed  that  the  observed  temperature 
may  be  somewhat  low  due  to  beat  losses  down  the  thermocouple  wire. 

This  possibility  will  be  investigated  experimentally  and  corrections 
applied  if  necessary. 
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Figure  3-  Typical  Theraocouple  Data 


R-2267-8 


18 


A  OlViaiON  oa  NOKTH  AMKNIOAN  AVIATIOKI,  INtt 


Two  conclusions  can  be  drawn  from  these  results: 

1.  The  uncertainty  in  the  x  measurement  prevents  the  differentia¬ 
tion  of  small  changes  in  surface  temperatures  as  would  be 
predicted  with  time,  axial  distance,  or  oxygen  flowrate.  Thus, 
much  of  the  discussion  on  this  point  in  the  previous  report 
becomes  academic. 

2.  All  of  the  extrapolated  surface  temperatures  are  extremely  low, 
unless  the  average  is  raised  well  over  200  C  by  a  heat  loss 
correction  which  appears  unlikely. 

The  isothermal  degradation  kinetics  for  FNHA  indicate  that  the  observed 
regression  rates  could  not  be  obtained  at  these  temperatures.  This  may  be 
an  indication  that  oxygen  or  other  chain- initiating  species  are  reaching 
the  surface  and  causing  wei^t  losses  far  in  excess  of  those  predicted 
by  pyrolysis  rate  data.  This  area  will  be  given  further  consideration. 
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FURTHER  THEORETICAL  ANALYSIS  OF  DIFFUSION 
MODEL  PROPOSED  PBEVIOUSLY 


The  experimental  study  of  hybrid  combustion  processes  has  concentrated 
on  the  burning  of  solid  fuel  with  gaseous  oxidizer  and,  specifically,  the 
burning  of  plexiglass  with  gaseous  oxygen.  The  experimental  results  have 
supported  the  idea  that  the  process  by  which  the  fuel  and  oxidizer  are 
brought  together  is  characterized  by  diffusion.  A  simplified  theory  was 
developed  following  the  Burke-Schumann  approach  (Ref.  6)  for  the  laminar 
diffusion  flame.  With  the  high  linear  flow  of  oxygen  in  the  restricted 
cylindrical  passage,  turbulence  must  dominate  the  pipe  flow;  and  therefore, 
it  was  necessary  that  the  parameters  of  the  theory  be  defined  for  turbu¬ 
lent  diffusion.  Thus,  no  distinction  is  made  between  turbulent  and  laminar 
flow  conditions  except  by  appropriately  chosen  parameters. 

Turbulence  and  diffusion  are  not  the  only  rate  controlling  factors.  The 
rate  of  combustion  is  limited  by  the  availability  of  fuel  through  the 
rate  of  gasification.  Turbulent  flow,  heat  transfer,  species  diffusion 
and  chemical  energy  release  all  play  coiig>licated  roles  in  hybrid  combus¬ 
tion.  In  an  effort  to  evaluate  these  effects  and,  possibly,  to  obtain  a 
simple  theory,  easy  to  understand  and  yet  generally  applicable,  it  seemed 
necessary  to  review  the  previous  analytical  work,  then  systematically  de¬ 
rive  and  reduce  the  differential  equations  involved,  in  the  hope  of  ob¬ 
taining  a  more  sophisticated  and  more  appropriate  model  than  the  one  pro¬ 
posed  \diich  was  obtained  throu|^  the  Burke-Schumann  approach. 
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The  first  model  was  obtained  previously  with  the  following  assumptions: 

1.  Chemical  reactions  are  confined  to  the  flame  front,  considered 
an  infinitely  thin  mathematical  surface  close  to  the  wall  to 
which  the  fuel  and  oxidizer  are  delivered  in  stoichiometric 
proportions  and  are,  therefore,  consumed  instantaneously.  For 
this  reason  no  fuel  exists  in  the  core  and  no  oxygen  exists  in 
the  annular  region  between  the  flame  front  and  the  wall. 

2.  The  pressure  gradient  is  negligible  and  the  linear  flow  vel¬ 
ocities  are  constant.  Hence,  the  axial  velocity  term  is 
neglected,  and  only  radial  transport  by  diffusion  is  considered. 
This  has  been  assumed  to  be,  primarily,  concentration 
diffusion. 

3.  Gas  mixtures  are  considered  to  be  binary  mixtures  composed  of 
the  major  constituent  (either  fuel  or  oxidizer,  depending  on 
location  beneath  or  above  flame  front)  and  the  rest  of  the  gas 
mixture  as  an  average. 

4.  The  product  of  gas  density  and  binary  diffusion  coefficient  is 
assumed  independent  of  radial  and  axial  distances  and 
temperature . 

5.  No  distinction  is  made  between  types  of  flow,  and  appropriate 
applicable  transport  coefficients  are  to  be  utilized  for  the 
axisymmetric  conditions. 

6.  The  binary  diffusion  coefficients  of  fuel  and  oxygen  through 
nonreacting  gas  are  numerically  equal. 
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Applying  these  assumptions,  the  continuity  equation  for  the  species  can 
be  simplified  and  combined  with  the  boundary  conditions  to  provide  a 
solvable  boundary-value  problem,  and  hence,  a  simplified  theory  of  hybrid 
combustion  as  summarized  in  Ref.  1. 

The  fundamental  equations  in  aerothermochemistry  as  proposed  by  von 
Karman  (Ref.  7)  can  be  the  foundation  for  a  more  sophisticated  theoretical 
development.  It  was  necessary  to  simplify  these  most  general  equations 
through  proper  consideration  of  the  conditions  of  this  specific  problem. 

The  following  assumptions  and  processes  (basically  the  same  as  those  in 
the  first  model)  were  then  applied  chronologically  to  obtain  a  solution: 

1.  The  compressible  gas  is  considered  to  be  in  thermodynamic  equi¬ 
librium,  i.e.,  the  concentration  of  each  species  is  a  unique 
function  of  the  local  pressure  and  temperature 

2.  Multicomponent  gas  mixtures  are  represented  by  the  fictitious 
binary  system  of  the  major  component  and  the  other  components 
by  average  characteristics.  The  specie  continuity  equation  is 
completely  described  by  one  equation  in  terms  of  the  mass  frac¬ 
tion  of  the  major  constituent  alone.  This  eliminates  chemical 
kinetic  considerations 

3.  All  parameters  are  independent  of  the  aximuthal  angle,  and 
azisymmetric  conditions  exist 

4.  Pressure  gradient  is  negligible 
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5.  By  analogy  with  the  Denison  (Ref.  8)  boundary  layer  equations 
for  chemical  ablating  surfaces,  the  applicable  pipe  flow  equa¬ 
tions  in  mean  quantities  can  be  written 

6.  The  cylinder  radius  is  observed  not  to  be  a  true  time  variable 
but  a  parameter  that  takes  on  definite  constant  values  at  each 
finite  burning  time 

7.  Existence  of  a  flame  front  \i4iich  is  the  site  of  all  chemical 
reaction,  i.e.,  an  aerotbermochemical  discontinuity 

Through  the  use  of  Eulerian  coordinates,  and  suitable  transformations  and 
decoupling  procedures  the  fuel,  oxygen,  and  temperature  profiles  can  he 
estimated.  Approximate  solutions  for  the  regression  rate  m  can  he  ob¬ 
tained  through  a  smthod  of  iteration  in  which  the  zeroth  order  corresponds 
to  the  first  approach  to  the  sisple  theory  developed  earlier.  Subsequent 
estimates  for  k  are  of  the  fora 


z'  being  the  function  of  the  axial  length.  This  should  yield  a  seasure  of 
the  modifying  effect  on  the  burning  rate  of  the  product  of  density  and 
axial  velocity  (£>u)  which  was  neglected  in  the  Burke -Schumann  iqiproach. 

Numerical  calculations  demonstrate  that  the  additive  term  in  the  above 
equation  decreases  the  k  toward  the  exit  in  contrast  to  the  experimental 
results.  It  is  believed  that  the  inaccuracies  in  the  theoretical  model 
are  inherent  in  the  assumption  Uiat  the  reactions  between  o:i^gen  and  fuel 
occur  only  at  the  infinitely  thin  flame  front,  and  that  the  reactants  are 
in  stoichiometric  proportions.  Physically  this  assumption  is  oversisplified 
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and  in  the  extreme  case,  after  the  oxygen  has  been  substantially  depleted, 
the  fuel  will  simply  gasify  because  of  the  hot  gas  near  it. 

Because  of  the  difficulties  in  defining  the  experimental  quantities  uti¬ 
lized  in  these  equations  and  in  solving  the  differential  equations  without 
employing  simplifying  assumptions  which  significantly  affect  the  solutions, 
no  further  work  is  planned  in  this  area. 
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SUMMARY 


The  r  vs  t  data  taken  at  two  oxygen  flowrates  Re2_  =  36,000  and  55|000 

°  c 

were  best  fitted  by  the  empirical  equation  r  >  A  +  Bt  %diere  c  »  0.73  for 
t  <60  and  c  :=  1.0  for  t  >  60.  Differentiation  of  this  equation  yields 
an  equation  from  which  the  instantaneous  rates  can  be  calculated.  These 
rates  have  been  tabulated  as  a  function  of  z  at  both  constant  burning 
time  and  constant  r.  The  ratios  of  the  regression  rates  at  the  two  flow- 
rates  have  a  dependence  on  z  vdiicb  cannot  be  accounted  for  by  changes  in 
turbulence. 

The  experimental  results  of  the  thermocouple  measurements  could  be 

fitted  by  the  psuedo-steady-state  heat  conduction  equation,  T  -  T  exp 

8 

(-bx/k);  an  average  t  has  been  calculated  from  b.  Two  conclusions  can 
be  drawn  from  the  results:  (l)  The  uncertainty  of  the  linear  measurement, 

X,  prevents  differentiation  of  small  changes  in  surface  temperatures  as 
would  be  predicted  with  time,  axial  distance  or  oxygen  flowrate,  and  (2) 
all  extrapolated  temperatures  are  extremely  low  to  account  for  the  observed 
regression  rates  from  isothermal  pyrolysis  rates  of  FMIA.  This  may  in¬ 
dicate  oxygen  reaching  the  surface  and  initiating  polymer  degradation. 

The  reliability  of  the  tesiperature  measurements  will  be  investigated 
further. 

In  an  effort  to  evaluate  and,  possibly,  to  improve  the  simple  theory  for 
hybrid  combustion  obtained  via  the  Burke-Schuman  approach,  a  theoretical 
development  was  made,  starting  with  the  von  Karman  fundamental  equations  for 
aerothermochemistry.  The  resultant  equations  should  have  shown  the  modi¬ 
fying  effect  of  including  the  pvt  term  neglected  previously.  Ihimerical 
calculations  indicate  that  the  additive  term  decreases  calculated 
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regression  rate  toward  the  exit,  which  is  in  contrast  to  the  experimental 
results.  It  is  believed  that  the  prediction  of  a  decrease  in  mass  regres¬ 
sion  rates  is  inherent  in  the  flame  front  assumption. 
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FUTURE  EFFORT 


Temperature  measurementa  within  the  solid  fuel  will  be  made  with  thermo¬ 
couples  of  smaller  diameter  wire.  These  will  be  used  to  determine  the 
error  caused  in  the  temperature  measurements  by  the  heat  losses  down  the 
thermocuple  wire  and  as  a  result  of  using  a  relatively,  large  bead. 

Burning  rate  data  at  increased  pressures  will  be  obtained.  To  check  the 
effect  of  radiation  absorption  of  tbe  surface,  additional  burning  rate 
data  will  be  secured  with  fuel  blocks  of  0.23-  and  0.30-percent  carbon 
in  PM4A. 
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